Quantum Depinning Transition of Quantum Hall Stripes 
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We examine the effect of disorder on the electromagnetic response of quantum Hall stripes using 
an effective elastic theory to describe their low-energy dynamics, and replicas and the Gaussian 
variational method to handle disorder effects. Within our model we demonstrate the existence 
of a depinning transition at a critical partial Landau level filling factor Av c . For Av < Av c , the 
pinned state is realized in a replica symmetry breaking (RSB) solution, and the frequency-dependent 
conductivities in both perpendicular and parallel to the stripes show resonant peaks. These peaks 
shift to zero frequency as Av — > Av c - For Av > Av c , we find a partial RSB (PRSB) solution 
in which there is free sliding only along the stripe direction. The transition is analogous to the 
Kosterlitz-Thouless phase transition. 

PACS numbers: 73.43.Nq, 73.43.Lp, 73.43.Qt 



I 

-a 
c 

o 
o 



> 

o 
in 
o 

m 
O 
-i— > 



I 

C 

O 

o 



X 



There is strong evidence from DC transport experi- 
ments [1] that states with stripe ordering form at cer- 
tain fillings of higher Landau levels (N > 2) [2]. These 
states exhibit a highly anisotropic, and apparently metal- 
lic, conductivity [1]. At low temperatures, as the partial 
filling of the highest occupied Landau level, Av, moves 
away from 1/2, the electrons in this level cease to con- 
tribute to the transport properties, and the system be- 
haves in a way typical of the integer quantum Hall effect 
[3]. One likely interpretation of this change in behavior 
is that the electrons in the partially filled Landau level 
become pinned by disorder when Av < Av c . The nature 
of this transition is the subject of this work. 

Microwave absorption measurements [4] provide addi- 
tional information about these systems. These experi- 
ments probe the dynamical conductivity of the system, 
a a f}(L>), which, in pinned systems, exhibits a peak at a 
frequency determined by the effective restoring force due 
to the disorder [5]. Existing data [4, 6] are suggestive of 
such a peak moving to zero frequency as the transition 
is approached from below, consistent with qualitative ex- 
pectations for a quantum depinning transition [7] . Inter- 
esting filling factor dependences of this peak have also 
been observed in Landau levels where a Wigner crystal 
is presumably pinned [8]. 

To examine the possibility of a depinning transition in 
the stripe state, we calculate the frequency-dependent 
conductivity using the replica trick and the Gaussian 
variational method (GVM), first introduced by Mezard 
and Parisi [9] and further developed by Giamarchi, Le 
Doussal, and their coworkers [10-13]. In the replica ap- 
proach, a pinned state is represented by one in which 
there is replica symmetry breaking (RSB). We demon- 
strate that for the stripe system, this leads to peaks both 
in Reo~ xx (u>) and in Hea yy ((v), albeit at slightly different 
frequencies and with different line shapes. (Here and in 
what follows, we adopt coordinates such that the stripes 
lie along the y direction.) An example of our results 



is illustrated in Fig. 1. A prominent feature of the re- 
sult is that the peak positions move to zero frequency as 
Av — > Av c from below. 

For Av > Av c , we find a different type of state in 
which the system is pinned for motion perpendicular to 
the stripes, but is free to slide along them. We call 
this a partial replica symmetry breaking (PRSB) state. 
The PRSB state has a number of striking properties, in- 
cluding a power law dependence of Kea xx (uj) ~ w 7 as a 
function of frequency with 7 continuously increasing with 
Av — Av c ; and a superconducting response [14] at zero 
frequency in Kea yy (uj), followed by an incoherent metal- 
lic response at finite frequency. The transition is in the 
Kosterlitz-Thouless universality class [7, 15]. It exhibits 
a jump in the low-frequency exponent 7 at the transi- 
tion, analogous of the universal jump in the stiffness of a 
thin-film superfluid and in the critical exponent of corre- 
lation functions [16]. The possibility of observing these 
properties in quantum Hall stripes and other analogous 
systems is discussed below. 

Model and method. We start with an action for an 
elastic system in a magnetic field [7] to describe the low- 
energy degrees of freedom of the quantum Hall stripes 
and their nonlinear coupling to the disorder, 
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Green's function of u in the pure limit (e xy — —e yx — 
l>£xx = e yy = 0), Ib is the magnetic length, and D a p(q) 
is the dynamical matrix. This last quantity is determined 
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[7] by matching the electron density-density correlation 
function obtained from the elastic model with that com- 
puted from microscopic time-dependent Hartree-Fock 
(HF) calculations [17]. In the low-energy sector, D a p(q) 
has a smectic form, D xx (q) ~ d xx (q x ) + n h q y , D xy (q) = 
D yx {q) ~ d xy (q x )q y , and D yy (q) ~ d yy {q x )ql. We note 
that alternative estimates of D were made by using an 
edge state model for the stripe system [18], which leads 
to different results than ours. We will comment on this 
difference below. The disorder potential V(r) in Eq. (3) 
is assumed to be Gaussian distributed with zero aver- 
age, V(r)V(r') = V 2 a x a y 5(r — r'), where a x and a y are 
the lattice constants of the stripe crystal. In Eq. (3) 
n(r,r) is the electron density operator whose Fourier 
transform we approximate by [10] n(q, r) ~ no{l — iq • 
u(q,r) +Ek/o/ dre iK -[ r - u ( r > r )]- ic i- r } where n„ is the 
average electron density and K is a stripe crystal recip- 
rocal lattice vector. For simplicity, we drop the iq-u term 
in n(q, t) since this cannot pin the electron system [10], 
and we keep only the lowest harmonics of the reciprocal 
lattice vectors, K x = 0,±2ir/a x and K y = 0, ±27r/a y . 
These approximations should not qualitatively change 
our results. 

We use replicas and the GVM to handle the dis- 
order effects. This method has been employed with 
great success in a number of systems [10, 11], although 
some controversy has arisen over previous applications 
to quantum Hall systems, for reasons that do not ap- 
ply to our present study [13, 19]. We introduce n 
replicas of the system where eventually we will take 
n — > 0. This involves setting S — ► Sq S ^ = E™=i 
with a the replica index; averaging over disorder intro- 
duces coupling among the replicas in the form S^J ~ 

-Vimp E",6=l Io /T dTl dT2 / dT E K #0 C0S{K- [U a (r, Tl)- 

u fc (r,T 2 )]}, with v imp = V 2 a 2 x a 2 y n 2 . The GVM [9] con- 
sists of replacing S^ 5 ' + S^J with a variational Gaus- 
sian action, characterized by variational Green's function 
G^(q, ioj n ) which is chosen to minimize the free energy. 
This leads to a set of saddle point equations (SPE's). 
It is convenient to parameterize the matrix G in terms 
of a self-energy matrix (, such that (G _1 )^(q, iu n ) — 

^aj -1 ^'^") S ab - (a b p{iuj n ). Note that we can safely 
assume £ has no q dependence because none emerges in 
the SPE's. We also set Ct = Ct = since it is a valid 
solution of the SPE's, and preserves reflection symmetry. 

Our SPE's are a natural generalization of those found 
using the replica and GVM on isotropic systems [10-12]. 
Here, we will present only the final equations, deferring 
details to a future publication [20] . After taking the limit 
n — > 0, the self-energy is characterized by a replica diago- 
nal component, Caa(^-Vi) — > Ca(i&n) and an off-diagonal 
function, (afi^n) — ► Ca( u )<L„,o, where < u < 1. It 
is the first of these that directly enters the frequency- 
dependent conductivity, while the second of these de- 
termines whether the system is in a replica symmetric 



state (constant ((u)), a RSB state (C(«) varying with 
u) or some other state. The dynamical conductivity is 
obtained [10] via 

a a p{uj) = {e 2 /a x a y )iwG T $(q = 0,w), (4) 

with G^(q,w) = G afi {q, iuj n -> u + iS), where G Q/3 
is the n — > limit of G^. The self energy Ca et ( w ) i s 
the analytic continuation of ( a {iu n ), which satisfies the 
following SPE's [20]: 

G?(w) = e a -[F a {u)-F a (0+)\, (5) 
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where e a = J du( a (u) — F a (0 + ) — 
2«im P Ek^o K 2 a J 1/T dr cxp{-T E M K E qiWB [1 
cosw„t] G w (q,w„)} and A^f) = E q ImG^(q, /). 

The constants e x , e y may be regarded as a measure of 
the strength of pinning by the disorder potential. When 
one or both of these vanish, the effective elastic constants 
G _1 (q = 0, uj n ) allow the system to be shifted as a whole 
without any energy cost. If one knows these constants 
then Eqs. (5) and (6) form a closed set of equations which 
may be solved numerically. A common further "semi- 
classical" approximation [10] involves expanding Eq. (6) 
for small A^, which is valid when the fluctuations of the 
stripes are small. While it greatly simplifies the numerics, 
this approximation is invalid near the depinning transi- 
tion since the fluctuations become arbitrarily large. We 
thus leave Eq. (6) in its present form. 

We are left with the task of computing e x and e y . As 
has been discussed elsewhere [10], this can be accom- 
plished without fully solving the SPE's by imposing the 
condition Re a ~ oj 2 for small u, which guarantees that 
the collective mode density of state vanishes at zero fre- 
quency. (This constraint can also be justified by requiring 
marginal stability of the replicon mode [10].) The second 
constraint can be found from the SPE's for ( a (u) with 
the assumption of a one-step RSB, which is a common 
solution in low-dimensional systems [10-12]. This leads 
after some work [20] to the condition 

ey/e x =Y, K >- Wm l[Y. K ^- Wm l ( 7 ) 

where W(K) = ± ^ K% / °° df A^(f) are Debye- Waller 
factors. These play a prominent role in the depinning 
transition: as we shall see, W(K) diverges whenever K 
has a component along the y axis as Av c is approached 
from below, leading to a suppression of e y . This behavior 
cannot be captured by the semiclassical approximation. 

Results - RSB solution. We have solved numerically 
Eqs. (5) and (6) along with the two constraints using an 



3 




0.025 



CO l B /e 



FIG. 1: Real parts of conductivities as functions of frequency 
in the pinned state (a) perpendicular to the stripes, (b) 
along the stripes. The U = 1 unit and «j mp = 0.0005e 4 /^s 
are used. In (a) all curves start from Kea xx = at 
u) — 0, and curves except for Av — 0.36 are lifted upward 
for clarity. Curves from right to left in (b) correspond to 
Av = 0.36,0.38,0.4,0.42,0.43,0.44,0.45,0.452,0.454, respec- 
tively. Inset in (b): peak positions Q px and fl py , in units of 
e 2 /Ib, as functions of Av. 



iterative method [20]. In what follows we present some 
results for electrons in the N = 3 Landau level, with 
a disorder level V[ mp — 0.0005e 4 /Z^. This is likely to be 
somewhat larger than experimental values, but we choose 
it for numerical convenience and do not expect our results 
to qualitatively change with smaller disorder strengths. 
The dynamical conductivities for these parameters when 
the system is in a pinned (RSB) phase are presented in 
Fig. 1. For Av well below Av c w 0.46, Rc<r xx has a pin- 
ning peak whose lineshapc is qualitatively similar to what 
is found using the semiclassical approximation [12]. The 
prominent behavior visible in Fig. 1 (a) is the collapse 
of the peak frequency Q px — > as the depinning tran- 
sition is approached. Experimental observations are so 
far consistent with this [4, 6]. Rea yy also has a collapsing 
peak, but the observed lineshapc is more interesting. Be- 
low the peak frequency in the range e y < u> < fl py 



the conductivity appears to tend toward a non- vanishing 
value when Av is sufficiently below Av c . The quantity e y 
turns out to be rather small due to a large Debye- Waller 
factor, and in this frequency range the system displays 
a behavior similar to the incoherent metal response at 
non- vanishing frequency of the depinned (PRSB) phase 
which we discuss below. For u -C e y , Rea yy (uj) vanishes 
quadratically with lo (not visible on the scale of Fig. 1) , as 
required for a pinned state. As Av — ► Av c , we eventually 
reach a situation in which e y and fl py are of similar order, 
in which case the pinning peak sharpens and grows quite 
large. This peak continuously evolves into a ^-function at 
zero frequency as the system enters into the PRSB state, 
so that the transition from pinned to depinned behavior 
is very continuous. 

Results - PRSB solution. For Av > Av c the state is 
characterized by e x ^ but e y = 0. This corresponds to 
a RSB solution for ( x (u) but a replica symmetric solu- 
tion for ( y (u). We call this the partial RSB (PRSB) state. 
In this situation, the system is pinned perpendicular to 
the stripe direction, but is free to slide along it. This 
is consistent with the results from a perturbative renor- 
malization group study of the same model [7] , where the 
coupling of the disorder to the motion of the system par- 
allel to the stripes was shown to be irrelevant when Av 
is close enough to 1/2. This irrelevance suggests that the 
PRSB phase should be in a superconducting state. This 
observation is born out by the presence of a <5-function 
m R,CtT yy £lt zero frequency. Remarkably, Reu^ vanishes 
at zero frequency, so that we find the PRSB state is one 
with an infinite anisotropy in the DC conductivity. This 
is not observed in DC transport experiments [1], and we 
comment below on what is missing from our model that 
leads to this discrepancy. The possibility of such behavior 
for quantum Hall stripes was first suggested in Ref. [14]. 

The origin of the PRSB and its structure may be un- 
derstood from the SPE (5). The state is characterized by 
e x ^ and lm( x ct (uj) ~ w at small cu, and CJ/ 6t ( w ) van_ 
ishing faster than linearly inw. It is easy to show in this 
situation, A y (f) ~ 1// at small /. After some algebra, 
we find [20] that when 
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2 > 1, (8) 



a self-consistent solution of Eq. 



(5) emerges, with 
u" and im (, y ~°(L0) ~ w 7+1 . The large 
value of this last exponent leads both to a power law 
dependence of Recr^ ~ w 7 and the 5-function response 



Re Q et (u) ~ to 2 and Im C y e *M 



It is also easy to show that 
so that, at low but finite 



at zero frequency in Recr OT . 
lim^^o Recr yy = const. > 
frequency the response is metallic. This combination of 
infinite and "incoherent" metallic response is typical of 
superconductors [21]. 

In the vanishing disorder limit, the minimum value 
of 7 from the formula (8) for which we can obtain a 
PRSB solution agrees precisely with the condition found 
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in Ref. 7 for pinning along the stripes to become irrele- 
vant. One remarkable consequence of this limiting value 
is that Rea xx ~ with 7 — > 1 as the transition is ap- 
proached from above, whereas just below the transition 
we expect, in the pinned state, Rc<r xx ~ uo 2 . Thus, the 
low frequency exponent jumps at the transition, in a way 
that is analogous to the universal jump in the superfluid 
stiffness and the critical exponent of correlation functions 
of the Kosterlitz-Thouless transition [16]. This behavior 
is similar to what happens in the roughening transition 
[15, 22]. 

In real DC transport experiments, one observes a fi- 
nite anisotropy rather than the infinite one found in the 
PRSB state. The missing ingredients from our model are 
processes allowing hopping of electrons between stripes. 
These processes are very difficult to incorporate into an 
elastic model. It is clear that, if relevant, such processes 
can broaden the <5-function response to yield anisotropic 
metallic behavior. Our results should apply at frequency 
scales above this broadening. Indeed, microwave absorp- 
tion experiments become quite challenging at low fre- 
quencies, and it is unclear whether existing measure- 
ments of the dynamical conductivity can access the low 
frequency conductivity in the unpinned state, whether 
or not it is broadened. In any case, it is interesting to 
speculate that a true 5-function response might be acces- 
sible in structured environments where barriers between 
stripes may suppress electron hopping among stripes [23] , 
or that there may be analogous states for layered 2+1 
dimensional classical systems of long string-like objects, 
which has been shown [15] to be closely related to the 
two-dimensional quantum stripe problem. 

The quantum depinning transition we find is unlikely 
to occur in models which preserve particle-hole symmetry 
(PHS) at Av = 1/2 [18]. Our model overcomes this 
limitation because the HF state we use spontaneously 
breaks PHS at this filling to arrive at a lower energy 
state than the simpler "box-filled" state [2] which has 
been used in the edge state description of the quantum 
Hall smcctic [18]. It is at present unclear if quantum 
fluctuations restore PHS to the quantum Hall smectic at 
Av = 1/2. Our results offer a falsifiable experimental 
test that can settle this question. 
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